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ABSTRACT 



^- . We present multi-frequency radio observations from the afterglow of 

/-^ . GRB 980519 beginning 7.2 hours after the gamma-ray burst and ending 63 

days later. The fast decline in the optical and X-ray light curves for this 



burst has been interpreted either as afterglow emission originating from a 
coUimated outflow - a jet - or the result of a blast wave propagating into a 



O ! medium whose density is shaped by the wind of an evolved massive star. These 

Q ! two models predict divergent behavior for the radio afterglow, and therefore, 



radio observations are capable, in principle, of discriminating between the two. 



0\ '. We show that a wind model describes the subsequent evolution of the radio 



afterglow rather well. However, we see strong modulation of the light curve, 
which we interpret as diffractive scintillation. These variations prevent us from 
decisively rejecting the jet model. 
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1. Introduction 



Both geometry and environment can affect the evolution of gamma-ray burst (GRB) 
afterglows ( |lVleszaros, Kees fc Wijers 1998| , [Panaitescu, Meszaros fc Kees 1998| ), giving us 
insight into such important issues as the total energetics, the burst event rate and the nature 
of GRB progenitors. The light curves and spectra from some of the earliest afterglows 
showed good agreement with fireball models of spherical blastwaves expanding into a 
homogeneous medium (e.g. Wijers, Rees fc Meszaros 1997 , Waxman 1997a ). However, 
over the past year, we have come to recognize a class of GRBs whose afterglows exhibited 
steeper than normal power- law decays (i.e. fi, oc t~^). Explanations for this behavior 
include invoking a large value for the energy spectral index {p ~ 4) of the radiation electrons 
( proot et ai 1998| ), a jet-like geometry for the relativistic shock ( pi.hoads 1997| , [Rhoads 



|1999| , |Sari, Piran fc Halpern 1999| ), and a class of inhomogeneous circumburst medium 
models, specifically those shaped by the winds (mass loss) of progenitor stars ( |Vietri l^Wl 
Chevalier fc Li 1999| , p~fc Chevalier l999| ). 



Empirical evidence connecting GRBs to the collapse of massive stars is accumulating 
from studies of GRB 980425 QGalama et ai 1998| , [Kulkarni et al. 1998| ), GRB 980326 
( IBloom et at. 1999b|) and GRB 970228 (|Reichart 1999| , [Galama et at. 1999| ). Further 
evidence comes from the studies of GRB environments, and includes the small measured 
GRB/host galaxy offsets ( [Bloom et ai 1999a| , [Bloom, Sigurdsson fc Pols 1999| ) and the 



growing number of optically obscured afterglows such as GRB 980329 ([Taylor et al. 1998|) 
and GRB 970828 (Djorgovski et al. 1999). 

The two models - the jet and the wind-shaped circumburst medium (WCM) - are well 
motivated by sound physical and empirical considerations. Jets have been invoked in part 
due to their universality, but also specifically to account for the large inferred isotropic 
energy release in some GRBs (e.g. GRB 990123; Kulkarni et al. 1999). The broad-band 
steepening of the optical and radio light curves seen in GRB 990510 is explained very 
well by the jet model ([Stanek et al. 1999| , [Harrison et al. 1999[) . Likewise, if GRBs are 
the end points of massive stars then it is inevitable that the circumburst medium will be 
inhomogeneous and reflect the mass loss history of the progenitor star. 

GRB 980519 occupies a central position in the debate of jets versus WCM models. It 
is the second brightest GRB in the BeppoSAX sample, and its X-ray and optical afterglow 
showed rapid fading, / ex t^^ (see pdalpern et al. 1999[ ). Based on the optical and X-ray 
data Sari, Piran fc Halpern (1999) advocated a jet model whereas Chevalier and Li (1999) 
advocate a WCM model. We discuss here the radio data in view of these two very different 
models. 
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2. Observations 



GRB 980519 was detected on 1998 May 19.51 UT with tlie Wide Field Camera (WFC) 
on board the BeppoSAX satelhte ( [Muller et at. 19^5] ) . A follow up observation with the 



Narrow Field Instruments (NFl) some 10 hrs after the burst detected a rapidly fading 
X-ray source ([Nicastro et al. 1998| ) with a power-law decay slope of — 1.83±0.3 ([Nicastro ei 



al. 1999| ). Within the WFC error circle, and coincident with the subsequent NFI location, 
Jaunsen et al. (1998) noted a fading optical counterpart. Djorgovski et al. (1998) showed 
that the optical transient (OT) exhibited a power-law decline with a ~ —2; here f{t) oc t". 
A summary of the X-ray and optical afterglow light curves can be found in Nicastro et al. 
(1999) and Halpern et al. (1999), respectively. 

We began observations of GRB 980519 with the Very Large Array (VLA) some 7.2 
hours after the initial 7-ray burst. The source was clearly detected on May 22; hereafter, we 
will refer to the radio afterglow as VLA J232221. 5+771543. A log of the VLA observations 
and a summary of the results can be found in Table |^. Except on the first day, when the 
bandwidth was halved in order to image the full WFC error circle, a 100 MHz bandwidth 
was used and all four Stokes parameters were measured. As these observations were begun 
when the VLA was in its most extended configuration and during the stormy summer 
season, considerable care was taken to accurately track the atmospheric phase variations 
across the array. Each scan on the source lasted three to five minutes and was bracketed on 
either side with a nearby phase calibrator (J0017+815). As a check on the integrity of the 
phase solutions a second phase calibrator (J2344+824) was observed frequently. The flux 
scale was tied to one of the sources J0137+331, J0542+498, or J1331+305. 

In addition to the VLA observations, we observed the OT with the Owens Valley 
Radio Observatory (OVRO) six-element array on the evening of 1998 May 21. GRB 980519 
was observed for 6.5 hr at a central frequency of 100 GHz with a 2 GHz bandwidth. The 
synthesized beam was 5"x4"; see Shepherd et al. (1998) for the procedure used to analyze 
the data. The resulting map had an rms noise level of 0.9 mJy beam~^ and the flux at the 
position of the OT was —1.2 mJy beam~^. Thus we place an upper limit (mean plus 3-cr) 
of 1.5 mJy beam~^. 

The light curve of VLA J232221. 5+771543 can be found in Figure |I|. The best fit 
position of VLA J232221. 5+771543, determined by a Gaussian fit of the combined 8.46 GHz 
data from May 22, May 24, June 2 and June 5 is (epoch J2000) a = 23''22'"21.50" (±0.02"), 
6 = +77°15'43.25" (±0.06"). The radio source coincides with the OT within the errors of 
the optical astrometry ( [Jaunsen et al. 199^ ). 
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3. The Light Curve: IntersteUar Scattering and Scintillation 

The light curve of VLAJ232221. 5+771543, with a gradual rise to a plateau followed 
by a decline below detectability around day 60, is qualitatively similar to that of many 
previously studied afterglows. Even though the data are sparse it is clear that there are 
some large amplitude variations (e.g. around about June 7), especially at 4.86 GHz. One 
could attribute the variations to residual uncorrected phases caused by poor weather. 
Fortunately, there exists a field source J232137.6+7715.0, some 2.5' from GRB 980519. As 
can been seen from the two lower panels in Fig. |l], the flux of this source is stable and 
does not show any suppression on those days when the afterglow was not detected. Phase 
incoherence is a multiplicative error and the stability of J232137.6+7715.0 is an empirical 
conflrmation of the robustness of our calibration procedure. Parenthetically we note that 
the fleld source is beyond the delay beam and is thus expected to suffer more photometric 
errors than a source close to the phase center. 

We believe that the large variations seen in the light curve of VLAJ232221. 5+771543 
are real, and following Goodman (1997) and Frail et al. (1997) we attribute these variations 
to interstellar scattering and scintillation (ISS). In particular, diffractive scintillation can 
induce extreme intensity variations (see below) and thereby account for the nulls in the light 
curve (Fig. [^). For example, in the relatively short time window of 12 to 19 days, we have 
four measurements with a mean of 103 /iJy and standard deviation of 103 /iJy at 8.46 GHz, 
i.e. 100% modulations. At 4.86 GHz we have mean of 80 fiJj and standard deviation of 
142 yLiJy, i.e. 180% modulations. The statistical fluctuations from thermal noise can explain 
about 30% and 40% of the flux modulations at 8.46 GHz and 4.86 GHz, respectively. For 
comparison, the fleld source had mean of 570 /iJy and standard deviation of 52 /xJy, while 
the statistical fluctuations are 85 /xJy. Therefore, the modulations in our data are likely 
dominated by scintillation rather than by statistical noise. We now proceed with the ISS 
hypothesis and use the observed variability to infer the angular size of the source. 

The strength of the scattering is reflected in the parameter SM, the so-called scattering 
measure. In the direction towards the GRB {l,b) = (117.96°, 15.26°), using the formulation 
of Taylor & Cordes (1993) we estimate SM = 5.8 x 10~^ m~^°/^ kpc; for reference, we note 
that this is a factor of two larger than towards GRB 970508 ([Frail et al. 19971) . We have 



assumed that the typical pathlength to the scattering medium in this direction is dscr=l-9 
kpc. Using this SM and Goodman's equations we expect to see strong scattering for 
frequencies below z/q = 15 GHz. The corresponding Fresnel angle and Fresnel size at z/q are, 
6fo = 1.5/iarcsec and Rfq = 4.2 x 10^° cm. 

For z/ > z/q, the scattering is weak and the modulations scale as (z//z/o)^^^/^^ < 1. The 
fact that we observe order of unity fluctuations at both 4.86 GHz and 8.46 GHz implies that 
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z/q > 8.5 GHz, compatible with the estimate from the SM. For u < uq we will see strong 
diffractive scintillation only if the size of the source, 9s < ^d = 9fo{^/^gT^^ ■ The other ISS 
parameters of interest are the decorrelation timescale (time for significant changes in the 
detected flux), tRiss = 3.9hr(z//z/o)^/^(i?Fo/4.2 x 10^°cm)/(t;x/30kms-^) and the bandwidth 
over which the diffractive ISS is decorrelated, Az/ = uq^u / uoY"^/^ . For z/q ~ 15 GHz the 
resulting {6D,tdm^ Az/) are (0.7/iarcsec, 1.8 hr, 900 MHz) at 8.46 GHz, and (0.4/xarcsec, 
0.9 hr, 80 MHz) at 4.86 GHz. In calculating tRiss we have assumed a transverse speed of 
30 km s~^ for the scintillation pattern across the line-of-sight. These estimates are rough 
in the sense that the estimate of SM (and hence z/q), for a given pathlength, can typically 
fluctuate by a factor of a few. Decreasing z/q has the effect of increasing Oo-, t^m and Az/. 

These strong diffractive fluctuations can be suppressed under two circumstances: (1) 
When the duration (At) and/or the bandwidth {B) of the observations exceed tRiss and Az/, 
respectively. The relevant figures of merit are given by rii = J B / ^v and n2 = -i/At/tRiss. 
(2) If the source size, Os-, exceeds Oo- The figure of merit here is ns = Os/Oo- If either 
Til or 71,2 > 1 then the observation will encompass more than one "scintel" (an island of 
constructive interference), thereby suppressing the variations. A large source will lower the 
modulation index to n^^ and increase the ISS timescale by the factor 713 ( [Narayan 1993|) . 



The fluctuations 5F of the observed flux F in the regime of diffractive scintillation will 
therefore be given by 

5F 

— = min(l, rii^) min(l, n^^) min(l, W^^). (1) 

F 

Our detection of strong diffractive scintillation at these frequencies suggests that all 
three figures ni,?T,2,?T-3 < 1; otherwise the fluctuations would be suppressed. From Table |1| 
and the estimate of the ISS parameters above we find ni(4.86 GHz) ~ 1, ?7,i(8.46 GHz) ~ 0.1 
and n2(4.86GHz) ~ 2, n2(8.46GHz) ~ 1. These estimates are compatible with the 
constraint discussed above. Interestingly, if z/Qwere any larger than 15 GHz, the fluctuations 
would have been suppressed since both rii and n2 would have become larger than unity at 
4.86 GHz. Our observations, therefore, limit vq to a narrow range of values (8 — 15 GHz) 
and likewise the estimates of ^z),^diflf, Az/. 

Thus, given the uncertainties involved, the strong modulations seen at both frequencies 
are quite reasonable. Furthermore, the strong modulations require that n^ be close to unity 
i.e. the source must be less than Od. Thus if our explanation that the strong variation seen 
in the light curve is due to diffractive ISS is correct, then the radio source must be less than 
0.4/iarcsec - a very small size indeed. Theoretical models predict a somewhat larger size. 
We will discuss this discrepancy in §|[ 

The strong modulations caused by scintillation do not allow us to accurately estimate 
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the spectral slope in the radio band. However, averaging over all observations from day 12 
on we find (3 = —0.45 ± 0.6 (where f^, oc z/^). This large range is compatible with both the 
early spectral index of 1/3 expected at low frequencies and with the late spectral index 
of ~ —1. It is inconsistent with the self-absorption spectral index of 2. The source was 
never detected at 1.43 GHz (see Table |l]). We conclude that VIA J232221. 5+771543 is an 
extremely compact (< 1/iarcsec) source with a self absorption frequency Uab between 1.43 
GHz and 4.86 GHz. 



4. Discussion 

The intense interest in GRB 980519 primarily stems from the fact that it was one of the 
first GRBs to be recognized to have a rapidly fading afterglow ([Halpern et al. 1999| ). Two 
entirely different models were proposed to account for the rapid fading: a jet expanding 
into a homogeneous medium ( ^ari, Piran fc Halpern 1999| ) and a spherical explosion into 
a wind-shaped circumburst medium (WCM; [Chevalier fc Li 1999| ). Both Sari et al. and 
Chevalier & Li used the same optical and X-ray data but ended up favoring two radically 
different models. 

The currently accepted view is that the afterglow emission arises in the forward shock 
of relativistically moving material ejected from a compact source (e.g. [Meszaros fc Rees 



19971 , [Waxman 1997b|) . It is assumed that in this shock the electrons are accelerated 
to a power- law distribution of energies (index p) above some minimum energy, 7,„meC^. 
Furthermore, it is assumed that there is a suitably strong magnetic field in the post-shock 
region. The energy density of the electrons and the magnetic field is assumed to scale 
linearly with the energy density of the protons. Gyration of the electrons in the magnetic 
field then results in a broad-band spectrum characterized by three frequencies, Vab^ ^m and 
Vc and the flux at frequency v^^ f^ ( |Sari, Firan fc Marayan 1998| ). The slope above u^ 
is given by f^ oc z/^^^^^)/^, and this steepens by one half at Uc due to radiative cooling. 
Below Um the spectral slope is the classical f^ oc v^^^, until it turns over below Uab due to 
synchrotron self- absorption, for which jy (xv^ . The dynamics of the expansion, the density 
distribution of the circumstellar matter and the details of the energy injection govern the 
the temporal evolution of these four parameters. However, the basic shape of the above 
spectrum does not change. Thus the key to distinguishing between competing models is 
through following the temporal evolution of the spectrum. 

Characterizing the broad-band spectrum by /(z^, t) oc f^v^ ^ we note from Halpern 
et al. (1999) in the first 1-2 days after the burst the following: Oopt = — 2.05±0.04, 
/?opt = -1.20 ±0.25, ax = -1.83 ±0.30 and /?opt-x = -1.05 ±0.10. Within errors the X-ray 
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afterglow declines as rapidly as the better measured optical afterglow. Likewise, within 
errors, the optical spectral index, /3opt is the same as the better measured optical-X-ray 
spectral index, /3opt-x- Hereafter we will refer to Oopt and /?opt-x by a and [5. 

Each model has a specific closure relation between a and (3. We now consider each 
model in turn. (1) A spherical explosion in a constant density medium would require 
H = or H = 1/2 depending on whether the range of observing frequencies is below 
or above Uc] here S{a,l3) = a — 3/2(3. The measured H = —0.48 ± 0.16 rules out this 
model (Halpern et al. 1999, Sari, Piran, & Halpern 1999). (2) A jet model would need 
J(a, (3) = a — 2(5 = —1, depending on whether the range of observing frequencies is below 
or above z/^. We find J = 0.05 ±0.21 which is seemingly more consistent with a jet expanding 
into a constant density medium with both the optical and X-ray frequencies below z/^. Sari 
et al. (1999) and Halpern et al. (1999) favor a jet model but with Uc between the optical 
and X-ray bands, in order to explain (3 with a standard value of p ~ 2.4. The apparent 
inconsistency with the value of J is attributed, in this interpretation, to the theoretically 
very long transition ( Panaitescu, Meszaros fc Rees 1998| , [Rhoads 1999| , [Sari, Piran fc| 



Halpern 1999| ) to the jet asymptotic power-law. The observed value of a may therefore 



be less than the theoretical one, if measured on a finite time interval. (3) An impulsive 
explosion in wind-shaped circumstellar model would need W{a,P) = 2a — 3(3 = —1, 1 
depending on whether the range of observing frequencies is below or above Vc- We find 
W = —0.95 ± 0.3 which is consistent with this model provided that both the optical and 
X-ray frequencies are below u^. Chevalier & Li (1999) justify the WCM model on this basis. 
However, they do note that their estimated Uc on day 1 (when a and (3 were measured) is 
between the optical and X-ray bands and argue that given the scanty measurements and 
rough theory the inconsistency is not particularly worrisome. 

Thus both the jet model and the WCM model can account for the observed optical 
and X-ray afterglow observations including the rapid decay. In the WCM model, the rapid 
decline is due to the gas density in the wind falling off as the inverse square of the distance 
between the shock boundary and the burst location, together with a larger value of the 
electron powerlaw index p ^ 3. In the jet model, the decrease is due to geometry. Due to 
relativistic "beaming" , only a small portion of the shocked gas is visible to the observer, a 
solid angle F"^ where F is the bulk Lorentz factor of the shocked gas. When F falls below 
9j^, the inverse of the opening angle of the jet, the observer will notice the finite edge of 
the emitting surface, thereby leading to a deficit of emission and thus a steeper than normal 
decline. Moreover, at about the same time, the jet will begin to expand sidewise, further 
enhancing the decline. 

We now discuss how radio observations can help resolve the ambiguity of the choice of 



models. For the sake of completeness, we also include a discussion of the simplest afterglow 
model (spherical explosion in a homogeneous medium). 

Spherical, homogeneous model. In this model, Uab oc t°, Um oc t^^/^ and fm = t^- As 
discussed in §^, Uab <8.46 GHz< z/^; this is certainly a secure statement two days after 
the burst. In this frequency range, we expect f{t) oc t^^^ until v^ reaches 8.46 GHz; we 
estimate (in the framework of this model) this will happen around 140 d. This estimate 
relies on a knowledge of the basic spectral shape of the afterglow ( [Sari, Piran fc JNarayan 



1998|) and the near simultaneous optical and radio data on 1998 May 28.59 UT to derive 
Um and fm at this time and evolve it forward until z/m=8.46 GHz. Thereafter the radio flux 
will decay in a manner similar to the optical and X-ray flux. As can be seen from Figure ^ 
the predictions are clearly inconsistent with the observations. Even the strong modulations 
produced by the ISS are unable to resolve this discrepancy (see below). Thus we can reject 
the spherical model in an entirely independent way from the closure method. 

Jet Model. We begin by noting that the steep optical decline had already started at the 
time of the first afterglow observation (ti) - the I-band observation by Jaunsen et al. (1998) 
about 7.5 hr after the burst. Thus the epoch at which the jet-like geometry of the emitting 
surface becomes apparent to the observer, tj, is less than ti. The epoch of the I-band 
observations is close to our first radio observation (see Table ||). In this "jet-dominated" 
regime the radio afterglow emission is /(i^afe < i^ < I'm) oc t^^^^ and /(z/ < Uab) oc t'^ (Sari 
et al. 1999). However, in §|| we argued that VLAJ232221. 5-1-771543 is an absorbed source 
with Uab between 1.43 and 4.86 GHz. The self- absorption frequency evolves slowly with 
time, Uab oc t~^^^. Thus we expect 8.46 GHz to be in the optically thin regime and therefore 
the flux at this frequency should decrease as oc t~^/^. A t~^/^ fit to the data, does not give 
a reasonable x^- However, the errors in the flux density (as quoted in Table |1|) may be 
underestimating the real uncertainties for a strongly scintillating source. If we fit the data, 
allowing for strong diffractive scintillation in the measurement uncertainties (i.e. a = S) 
then a reasonable solution can be found S=210t^ yuJy {x^=Q with 13 dof), where td is 
time in days. The sharp break in the light curve at tm = 26 days corresponds to 1/^=8.46 
GHz and was derived by evolving the spectrum forward from t = 1.08 days (with /3=l/3 for 
u < Um and P = /5opt-x for u > Um), when high quality optical and radio data existed. 

Wind-shaped Circumstellar Medium Model. Using the contemporaneous X-ray, 
optical and radio data taken up to three days after the burst. Chevalier & Li (1999) fit the 
circumstellar model to derive the fireball parameters (total energy, energy power law index, 
fraction of energy in electrons and magnetic fields) and the wind properties, for an assumed 
redshift of z=l. In their derivation, the 8.46 GHz emission initially originates below Uab and 
Um but the subsequent evolution moves them into this band as Uab oc t~^^^, Um oc t~^/^ and 
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fm oc t~^/^. Thus in the circumstellar model one predicts an initial rise of the radio flux as 
the source opacity decreases, followed by a plateau when Uat <8.46 GHz< z/^ and then a 
steep decay when Um <8.46 GHz. We have plotted the model fit of Chevalier & Li (1999) 
(labeled "Wind") in Fig. 0. We stress that the model was derived on the first three days of 
afterglow data, spanning 8 orders of magnitude in frequency, and yet without any further 
adjustments their model fits the next 60 days of VLA measurements remarkably well. 

As we have noted in the previous section, the angular size of the afterglow is less 
than l/xarcsec even after about 15 days. In the three models, spherical, wind and jet, the 
afterglow angular size at that time is given by 

Osphere = 2.8 /xarcsec (^^) Djl,sEli\-'/\t/15 days)^/^ (2) 

ew^nd = 2.2 /iarcsec [-^) D-^%^Eg^A-^'\t/lb days)^/^ (3) 

Oj = 1.7/iarcsec (^) ' ' D^^,Eli\;'^'{t j/ Shi)' /'{t/ 15 days)'/'. (4) 

Here, Da,28 is the distance in units of 10^*^ cm, £'52 is the inferred "isotropic" energy released 
in the explosion in units of lO^^ergs, rii is the surrounding density in the constant density 
case. A^ characterizes the wind density as p{R) = 5 x 10^^ A^R~^ g cm~^, with as a Rem 
the wind radius in cm. These sizes are quite similar to each other. This stems from the 
very low dependence of the size as a function of density and from the fact that the density 
in front of the shock, at a time of 15 days, in a wind model with typical parameters is not 
very different from regular ISM densities. It can be seen that the small inferred size from 
the scintillation is challenging for all these models. It requires either a dense wind, or a jet 
with very small opening angle, that begins to spread at a very early time. It is hard to 
reduce the size below 1 /xarcsec for any reasonable choice of the parameters in the spherical 
constant density model. The agreement with scintillation theory would improve if either 
SM or the effective distance of the scattering screen dsa were lower by a factor of two. 
Along specific lines of sight the uncertainty in either of these quantities could vary by this 
amount due to the non-uniform nature of the ionized medium. In such a case the inferred 
size of 0.4 /iarcsec would increase by a factor of 2.5, enabling both jet and wind models to 
fit the data. 



5. Conclusions 

The observed steep decline in the optical and X-ray light curves of the afterglow from 
GRB 980519 lends itself to two equally compelling hypotheses, the first being that the 
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afterglow emission originates from a collimated outflow (or jet), and the second, that the 
emission is the result of a blast wave propagating into a medium whose density is shaped by 
the wind of an evolved massive star. These competing models predict divergent behavior 
for the evolution of the radio emission. In the jet model the radio flux is expected to 
decrease continuously after the jet edge becomes visible. In contrast, the circumstellar 
model, specifically the Wolf-Rayet wind model of Chevalier & Li (1999), predicts a linear 
rise, followed by a broad plateau and a late decay in the radio hght curve. 

In this case, scintillation, low-signal-to-noise and sparse data at early times, has limited 
us from decisively choosing either model. The influence of ISS could be reduced, in future 
afterglow observations, if one has more measurements at early times (the first two days is 
where the models differ the most). Such observations, separated by a few hours (a time 
which is larger than tjss), could be averaged to produce a more solid estimate of the emitted 
flux, free of large modulations. Nevertheless, this result illustrates the unique diagnostics 
provided by radio observations and its potential for unraveling the origin of GRBs. Previous 
examples include the demonstration of jets in GRB 990510 ( [Harrison et al. 1999|) and the 



discovery of a reverse shock in GRB 990123 ( Sari fc Piran 1999 , Kulkarni et at. 1999|) 
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Table 1. VLA Observations of GRB 980519 



Epoch 


At 


TOS 


S8.46 ± 0" 


S4.86 ± cr 


S field ± O- 


S1.43 ± f 


(UT) 


(days) 


(min) 


(/iJy) 


(MJy) 


(MJy) 


(MJy) 


1998 May 19.81 


0.30 


99 


49±28 








1998 May 20.59 


1.08 


54 


64±27 








1998 May 22.35 


2.84 


112 


103±19 








1998 May 24.96 


5.45 


112 


127±20 








1998 May 31.54 


12.03 


132 


40±25 


25±27 


591±80 




1998 Jun. 02.56 


14.05 


125 


142±29 


292±41 


607±100 




1998 Jun. 05.14 


16.63 


106 


230±31 


16±33 


596±80 




1998 Jun. 07.43 


18.92 


183 


1.2±32 


-13±37 


493±82 


-4.5±34 


1998 Jun. 11.94 


23.43 


69 


82±40 


215±44 


537±90 




1998 Jun. 18.58 


30.07 


78 


66±23 


5±31 


676±70 




1998 Jul. 07.26 


48.75 


100 


78±27 


133±42 


480±71 




1998 Jul. 20.23 


61.72 


124 


0.9±30 


190±40 


356±73 


15±42 


1998 Jul. 21.25 


62.74 


94 


42±32 


57±45 


480±79 


-35±58 



Note. — The columns are (left to right), (1) UT date of the start of each observation, (2) time elapsed 
since the GRB 980519 event, (3) total duration of the observing run, (4) 8.46 GHz flux density of the radio 
transient (RT), with the error given as the root mean square flux density, (5) 4.86 GHz flux density of the 
RT, (6) The integrated 4.86 GHz flux density of J232137.6+7715.0, a field source 2.5' from the GRB, and 
(7) 1.43 GHz flux density of the RT. 
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Fig. 1. — Observed radio lightcurves at 8.46 GHz and 4.86 GHz. The two top panels are 
light curves for the radio afterglow from GRB 980519, while the bottom panel is the light 
curve for J232137. 6+7715.0, a field source 2.5' from the GRB. Flux densities (and errors) 
are taken from Table |l|. 
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Fig. 2. — Observed and model lightcurves at 8.46 GHz. Detections are indicated by the filled 
squares. Upper limits for the non-detections (open triangles), are plotted as the measured 
flux at the position of the radio transient plus twice the rms noise. Three different model 
predictions are plotted. The circumstellar model "Wind" is taken directly from Chevalier & 
Li (1999) with no modifications. The basic adiabatic, forward shock model with spherical 
expansion into a homogeneous medium is indicated as "Homogeneous". The predicted 
behavior for the jet model of Sari et al. (1999) is given by "Jet". See text for additional 
details. 



